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Abstract 

One way to treat the infrared divergences of the electroweak Next-to-Leading-Order (NLO) dif- 
ferential cross sections to parity-violating (PV) electron-proton scattering is by adding soft-photon 
emission contribution. Although more physical, the results are left with a logarithmic dependence 
on the photon detector acceptance, which can only be eliminated by considering Hard Photon 
Bremsstrahlung (HPB) contribution. Here we present a treatment of HPB for PV electron-proton 
scattering. HPB differential cross sections for electron-proton scattering have been computed using 
the experimental values of nucleon form factors. The final results are expressed through kinematic 
parameters, making it possible to apply the computed PV HPB differential cross sections for the 
analysis of data of a range of current and proposed experiments. 
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I. INTRODUCTION 



Electroweak properties of the nucleon can be studied by parity- violating electron-nucleon 
scattering at low to medium energies Such experiments can measure the asymmetry 
factor coming from the difference between cross sections of left- and right-handed electrons. 
This asymmetry between left- and right-handed particles, as a result of a parity-violating 
interference between the weak and electromagnetic forces, is clearly predicted in the Standard 
Model of Particle Physics. 

Extracting the physics of interest from the measured asymmetry requires evaluating NLO 
contribution to electroweak scattering at very high precision. The method for evaluation 
of the electron-nucleon up to NLO differential cross sections most commonly found in the 
literature is to follow the Feynman rules for the particles of the Standard Model. The 
dominant contribution normally comes from the leading order (LO) correction in pertur- 
bation theory. Some of the electroweak NLO contributions to intermediate energy, parity 
non-conserving semi-leptonic neutral current interactions have been addressed previously in 
Ref. jfl] also estimated effects due to an intrinsic weak interaction in the nucleon 
(e. g. the anapole moment) in chiral perturbation theory, and found the anapole moment 
contributions insignificant, only slightly enhancing the axial vector NLO contribution. 

Later work in [7j took the advantage of the modern computational opportunities and im- 
proved the techniques for one-quark NLO computation by retaining analytical momentum- 
dependent expressions, and providing the numerical evaluations of 446 one-loop diagrams. 
It also included calculation of the soft photon emission contribution. However, in [7|, even 
after removing infra-red (IR) divergences through soft-photon emission corrections, calcu- 
lated one-quark NLO contribution show a logarithmic dependence on the detector's photon 
acceptance parameter AE. 

The article presented here demonstrates that elimination of this dependence can be 
achieved by adding the Hard-Photon Bremsstrahlung (HPB) differential cross section. We 
express general electroweak couplings by inserting appropriate form factors into vertices and 
construct HPB differential cross sections as a function of Mandelstam invariants. For each 
set of experimental constraints, integration over the emitted photon phase space can be 
performed numericall y. A nalytical results of this article can be used for several recent PV 
experiments @, |, MM Q 

The article provides a detailed description of both hard- and soft-photon emission treat- 
ment of infrared divergences the PV electroweak interference and pure weak contributions 
to the total differential cross section. As an example, we choose to consider electron-proton 
scattering, as one of the most relevant cases from a physics perspective. However, the same 
technique of treating infrared (IR) divergences can be expanded to neutron or any other 
baryon target if same effective structure of the coupling is used. 



II. SOFT-PHOTON BREMSSTRAHLUNG 



If the structure of the nucleon is investigated using weak, neutral current probe 0, II 
it is necessary to enhance the weak contribution in electron-nucleon scattering by exploiting 
the parity- violating nature of the weak interactions and constructing the following quantity 
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Here, A LO is a leading order asymmetry measured in the parts per million (ppm) and second 
term of Eq.(QQ) is a parity- violating NLO contribution to the asymmetry. 

Finally, if the axial- vector (Ci) or vector- axial (C2) form factors of the parity violating 
amplitude are studied [13], which can be taken from PV Hamiltonian 



h pv _ (f_F_ 



d (u e Y7 5 Ue) (uN^Un) + C 2 (Ue^Ue) (u N ^ & U N ) ] , (2) 



V2 

with perturbative expansion resulting in 

C li2 = Cg + C»£° + O (a 3 ) . (3) 



All of the above NLO contributions to the either asymmetry (Eq.(pQ)) or PV form-factor 



(Eq.jH])) in general can be infra-red divergent |l4j, and can be treated by the soft and 



hard-photon emission contribution shown on Fig. (DQ). 




Figure 1: Hard Photon Bremsstrahlung diagrams in electron-proton scattering. 
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The differential cross section associated with bremsstrahlung emission in electron-nucleon 
scattering can be described by the following formula, 
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The first term of Eq.((4j) is responsible for the cancellation of IR divergences if a parity 
conserving electromagnetic probe is used. The second term, when used in the asymmetry 
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is responsible for canceling IR divergences in Eq.(QQ). Finally, IR divergences in the NLO 
form- factors (Eq. ((3j)) are indirectly treated by the third term of Eq.Q, and will be discussed 
later in this article. 

Generally, bremsstrahlung diagrams can be described as 2 — > 3 processes in which inte- 
gration over emitted photon's phase space should be performed. If the momentum of the 
emitted photon is small enough to be neglected in the numerator algebra, we can present the 
bremsstrahlung cross section as a soft photon factor multiplied by the tree level differential 
cross section of 2 — ► 2 process. 

Let us consider an example. The scattering amplitude for the first diagram of Fig.(pp), 
for the neutral current reaction, has the following structure: 
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where the photon polarization vector enters as e a {k^)] T z _ e , ± Z _ N aim j- 7 _at 



r^-jv an d ^-n are the cou- 
plings of electron with Z boson, nucleon with Z boson, and photon with nucleon, respectively, 
defined as 
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Here q corresponds to the momentum transferred to the nucleon from the vector boson. The 
shortened notation of s\y and cw refers to sin 6^ and cos6 w , the Weinberg mixing angle. 
For the form factors fi^(q), ^1,2(9), and gi(q) we have used 
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Fi, 2 (q) = F li2 (0)G N (q) , 

(7) 

F$ n) (0) = (0) - Ff >2 (0) , < 2 (p) (0) = f* 2 (0) - Fl 2 (0) , 

9M = -^-G N (q), 

with Fi 7 2 (0) and ^ (0) defined as the nucleon's Dirac, Pauli, and axial form factors at zero 
momentum transfer and corresponds to the electric charge, anomalous magnetic moment, 
and axial charge, respectively. Here, we use the universal formfactor of monopole or dipole 
type G N (q) = (jtz?)" with n = 1, 2 and A 2 = 0.83m^. 

In the soft-photon emission limit the coupling between the emitted photon and the nu- 
cleon, r"_ 7V (fc 5 ), is just equal to ieQj a , where charge Q — 1(0) for proton (neutron). The 
numerator of nucleon's propagator /k 2 — /k 5 + m N can be replaced by /k 2 + m N , and 
{k 2 — ks) 2 — m 2 N can be easily simplified into —2 {k 2 • k§). Using the Dirac equation for free 
spinors, we have 
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= -leQ—r. — — w(mjv, fc 2 )) . 

Now we can present the soft-photon amplitude in the following form: 
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Here, M LO is a tree level amplitude of 2 — > 2 process. Eq.((9|) can be used for the other 
photon emission diagrams with a different factor n(ki, k 5 ) = ±ieQi ^..j^ ■ Now we can 
sum over all four graphs on Fig. (Op) and square the total amplitude to get the following: 
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The photon couples to a current which is conserved: k^M^ = 0. This fact, and the summa- 
tion over all photon polarizations gives us the possibility to replace ^e{kiY{kj) l/ e*^e v with 
—gfj,v{.kiY{kj) u = —(ki-kj). The last step is to integrate over the emitted photon phase space 
^ fc 5 = (27r3)2fc an d regularize the infrared divergence by assigning to the photon small rest 
mass A. This dependence on the rest mass of the photon will be canceled when added to IR 
divergent NLO contribution (see Eqs.(ED, Ej)). The resulting soft-photon emission differential 
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Table I: Soft photon emission integral parameters of Eq. (2.10) 
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Here AE is the maximum possible energy of emitted photon for which the soft-photon 
approximation is still valid. Numerical analysis leads to a typical constraint A < AE < 
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v y — 2a e -e ' an( l l^lj are ^ ne f erm i° n ' s energy and spatial momentum in center 
of mass reference frame, correspondingly. The parameter can be extracted from Table 

The dependence on AE is canceled as a result of adding soft and hard-photon emission 
differential cross sections. It is worthwhile to mention here that da so f t is proportional to 
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III. HARD-PHOTON BREMSSTRAHLUNG 



This section gives details on the evaluation of hard-photon bremsstrahlung differential 
cross section. The results are expressed in a form convenient for further analysis. 



A. Electron-Nucleon Scattering 



In the case where the energy of the emitted photon (k® > AE) can no longer be neglected 
in the numerator algebra, we have to account for all the differences between hard- and soft- 
photon emission. Besides the fact that the hard-photon amplitude will have k$ in the 
numerator, calculations for differential cross section will have to include helicity matrix 
elements with extended set of Mandelstam variables. These matrix elements come from the 
use of the momentum conservation law for 2^3 process. Thus, the helicity matrix elements 
will depend on the extended set of Mandelstam variables: 
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The evaluation of the interference term Re(M 2 ot 



Ml7z) and 



Ml7% 



is somewhat cum- 
bersome because it includes calculations of 3136 helicity matrix elements. Here, t' — m 2 z 
can be replaced by t — m| due to the fact that {t, t'} <C m 2 z . The evaluation of the HPB 
contribution can be further simplified by spliting the amplitude into two parts: 
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Here, M 2 ^ 7 j is the total amplitude with the momentum of the emitted photon k 5 removed 
from the numerator of Eq. (TT4l) . and M 2 ^ 7 } is everything that is left up to order 0(k 5 ). 
Now, the squared amplitude for the neutral current reaction has a simple form: 
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The interference term can be calculated as 
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Applying Dirac equation in Ml? z ^ amplitude, we can simplify our calculations consid- 
erably. The first terms of Eq. (Tl6l) and (fT7l ) can be obtained from 
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where 
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The term {f^z^-n) represents a coupling which was modified in a way so it would no 
longer have dependency on the hadron formfactor Gn (t), and no longer contain momentum 
of the photon in its Pauli part of coupling: 
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In Eq.(fT9l) and Eq.(fT8i). the coupling I^r ^ again was replaced by (ieQ) 7". It is straight- 
forward to see what after the integration over the phase space of the emitted photon only 
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will have a logarithmic dependence on the photon de- 
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and Re(M 2 



when combined 



with the soft-photon bremsstrahlung differential cross section, will be responsible for the 
cancellation of Log(AAE 2 ) term in Eqs.f fTTT) and (fT2l ). 

Further numerical analysis shows that, when integrated, the second and third terms of 
Eq.jHJ) and Eq.(TT7l) have no logarithmic dependence on AE. They both are small compared 
to the first term when energy of incident electrons is in the domain of the current or proposed 
PV experiments. This simplifies calculations of PV HPB contribution considerably, since 
the only first term of the Eq. (fT6l) and Eq. (fT8l) has to be considered in the calculations. Here 
we provide details on how to calculate first terms of Eq. (fT6l) and (PT7l) explicitly. Although 
details of calculations for the rest of the terms are not shown in this article we have them 
included in our numerical analysis. 
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The scalar products (ki ■ fa) are Lorentz invariants, and can be replaced with the Mandelstam 
variables as 
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As for 
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and Re(M' a M' a z ) l,r, we have detailed expressions given in the appendix 



of this article. The helicity matrix elements were computed with the help of FormCalc [15 



Now we are ready to proceed to the next sections, where we shall give the details on 
the parametrization of the emitted photon's phase space, and numerical details on the 
calculations of the PV HPB contribution. 
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IV. HPB DIFFERENTIAL CROSS SECTION 

Parametrization of the phase space for (2 — > 3) process has been chosen according to Fig. 




Figure 2: Phase space for the emitted hard photon. 

Here, the angle 9 is a scattering angle and £ corresponds to the angle between emitted 
photon and scattered electron. The momenta are represented as 



ki = {Ei, 0,0, Pi n }, 
k 2 = {E 2 , 0, 0, -p in } 
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where unit vectors are 
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e 5 = 



/ cos(6 l ) cos (77) sin(£) + sin(0) cos (£) 

sin(?7) sin (£) 
^cos(^) cos(£) — sin(6>) cos (77) sin(£) 



For on-shell particles, the incident momentum pi n can be found as 



with 
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The center-of-mass energy E cms can be determined as follows: 



E cms = \Jm 2 e + m 2 N + 2Ei ab m N . (29) 
Momentum k± is determined by the four-momentum conservation law in the cms frame: 
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and 



k 3 + k 4 + k 5 = 0. 
The HPB differential cross section reads as follows 
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where $ is a flux factor and given by 
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The (2 — ■> 3) process phase-space element c£T 3 is 
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mentum k 4 , we arrive at 
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with = d cos 9 d(f) and dQ$ = dcos£ drj. The remaining delta function 5 (k 2 — m 2 N ) will 
be used to eliminate integration over the scattered electron energy k®. 
We need to do some modifications first: 



Now, using 
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The electron mass can be considered as a small parameter with respect to k®. In this case, 



we replace 



by k®, so that 
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Substitution of Eq.(l39l) into Eq.(l38l leads to the following: 
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The property of the delta function 5 [g (k%)] = J2i | g ' 3 ( r -)| ( r * * s *~^h r00 ^ °^ the e Q ua tion 
k\ — m 2 N = 0, solved with respect to k®) makes it possible to replace 5 (k 2 — m%) by 
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The delta function S (k% — r) will eliminate integration over k® leaving k® = r. Integration 
over the emitted photon's phase space dk^dfl^ can be performed numerically using the cuts 
on the photon's energy k® : 



Finally, Eq. (|32|) becomes 
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leaving the final differential cross section differential with respect to the scattered electron 
solid angle dQ 3 . 
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A. Numerical Test 



Now we can introduce details of the treatment of infrared divergences in the parity- 
violating formfactors Ci® +NLO . Since the PV amplitude derived from the Hamiltonian 
Eq.(ED can be used in the calculations of the asymmetry Eq.(QQ), the differential cross section 
for the neutral current reaction can be computed as: 
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The contribution of the soft- and hard-photon bremsstrahlung modifies differential cross 
sections in the Eq.jl]) according to the following: 

d °R, L = d(7 R,L + { da L0^ 2 ) RL ■ <Wt + da z R f PB . (45) 

In order to combine HPB differential cross section with the soft-photon emission contribution 
factor 5 so ff, and all this with parity- violating formfactors Ci ° +NLO ', we propose the following 
parametrization for the HPB differential cross section: 

j Z,HPB i Z,HPB 

i~Z,HPB ( r Z,2^2\ aa R — acr L (i Z,2^2\ r ( AaX 

As can be easily seen, the substitution of Eq. fl46l) into the expression for asymmetry 
Eq.(DQ) will leave terms related to the neutral current reaction HPB in the usual form 
(da R HPB — daf i ' HPB ^j . It is worth noting that the term (da^ + da^) has a dominant con- 
tribution from the parity-conserving part of the differential cross section. Because of that, 
denominator of Eq. (pQ) is left without parity- violating soft- and hard-photon bremsstrahlung 
terms. Combining the soft term of Eq.JTT]) and the HPB term of Eq. (l46|) with PV formfac- 
tors, we can write 

riLO 

C™> = + Y (<W + Shpb) ■ (47) 

For the case of {e — N} scattering, we will show numerical contribution from the SPB and 
HPB terms, taking into account only the IR finite part of the soft-photon bremsstrahlung 
only. We can do so because IR divergences are canceled when {e — N} PV formfactors 
Ci2° are combined with the second part of the Eq. (1471) . Moreover we will treat formfactor 
Gn (?) using monopole approximation {n = 1) in our numerical tests. 

Let us start with demonstration that, indeed, we do not have a AE dependence in the 
term | (5 softer Shpb) for a kinematic point relevant to the Q wea k experiment. We take 
Ei ab = 1.165 GeV and Q 2 = 0.03 GeV 2 . During the numerical integration, we have used the 



adaptive Genz-Malik algorithm which is implemented in the Mathematica program 16j. For 
electron-proton scattering, the term | ($ so ft + $hpb^) for different values of AE is shown in 

the Table (JTTJ) . We see that the variation of | (S so f t + $hpb) is of order of ~ 0.1% which 
is coming from the statistical error of integration. The same can be done in the analysis 
of the AE dependence for the PV asymmetry due to soft and hard photon bremsstrahlung 
(see Eqs.((5!), (TTT1) and f|43l ). Data for table (IIH|) have been computed using the same 
integration technique, and it is clear that variations of the asymmetry are < 0.05 %. In the 
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AE {y/s) 



hoft 



+ &HPB 



10" 


-3 


-0.16966 


10" 


-4 


-0.16980 


10- 


-5 


-0.16982 


io- 


-6 


-0.16984 


10" 


-7 


-0.16982 



Table II: Dependence on the photon detector acceptance (electron-nucleon scattering case Ei a b 
1.165 GeV, Q 2 = 0.03 GeV 2 ) 



AE (y/s) A B (W- 8 /dT 2 ) 



io- 


-3 


1.00019 


io- 


-4 


1.00061 


io- 


-5 


1.00068 


io- 


-6 


1.00068 


10" 


-7 


1.00067 



Table III: Dependence of the bremsstrahlung asymmetry given in the units of 1/dT 2 = 4 ^^ | V ^ on 



the photon detector acceptance AE (electron-nucleon scattering case E iab = 3.0 GeV, Q 2 = 0.3 
GeV 2 ) 



test of independence from the AE parameter we took one of the kinematic points of the 
GO experiment, with Q 2 = 0.3 GeV 2 . For the complete analysis of e — p PV scattering 
asymmetries it is required to include all the LO and NLO contributions, which will be left 
to a future publication using the treatment of IR divergences described in this article. 



V. CONCLUSION 

The calculation routines, tested for electron-proton scattering and presented in our pre- 
vious work are valid for any electroweak processes involving particles of the Standard 
Model. The HPB contribution computed in the current work can be applied for virtually any 
scattering process. Again, when computing hard-photon bremsstrahlung terms for electron- 
proton scattering, we use effective Z — p and 7 — p couplings with monopole type form 
factors. 

The enormous size of the complete analytical expressions involved makes it impossible to 
present them in this paper. The complete analytical expression in the Mathematica file is 
available from authors upon request. 

We observed that for the energy range employed by the PV experiments the HPB differ- 
ential cross section is dominated by part of HPB amplitude without the photon momentum 
in the numerator. It is still necessary to keep in mind that the process is (2 — > 3) when the 
cross section is calculated. Terms proportional to O (k 5 ) tend to be important for higher en- 
ergies. This simplifies calculations of the HPB cross section for the considered experiments 
significantly, as it simplifies the numerator algebra. 

We split the amplitude in two parts, with one part being the amplitude without the mo- 



ll 



mentum of the emitted photon in the numerator. This step is important, because, according 
to the numerical analysis performed, this term has a strong dominant structure similar to 
the soft-photon emission factor. Another interesting result of this work is that all of the 
effects, including soft- and hard-photon bremsstrahlung terms, can be now accounted for on 
the level of PV formfactors. 

The proper account of the soft- and hard-photon bremsstrahlung effects has allowed us 
to achieve final results that are free from a logarithmic dependence on the detector photon 
acceptance parameter. 
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VI. APPENDIX 



Here we give the detailed 



Kz 



L,R 



expressions used for the calculations of 



Ml? 



L.R 



m 



Eq. (T22l for left-handed incident electrons: 



Kz 



2 _ 4a 3 7r 3 (1 - 2s 2 w f 

L ~ r 2 s 2 



(/ 2 (o)r 

Am 2 N 



'lQm% — 8m%(s + s' + u + u) + Am 2 N (s + u)(s + u') + 



t'(us - ss + tt' + su' - uu')) + g z R N f 2 (0)(4m% - Am 2 N (u + u ) - ss + tt' + s'u + sv! + 

3uu') + 2 (gt N ) 2 (ml - s)(m 2 N -s') + 2 (g z R N ) 2 (m 2 N - u)(m 2 N - u') + 

g z L - N (Ag z R - N m 2 N t + f 2 (0)(4m% - Am 2 N (s + s) + 3ss' + tt' + s'u + su - uu')) 
and for the right-handed incident electrons: 



(48) 



M' z 



16a 3 7r 3 s 2 



w 



R 



Am 2 N 



[16m% — 8m%(s + s' + u + u) + Am 2 N (s + u)(s' + u') + 



t'(us — ss' + tt' + su — uu')) + g R N f2(0)(Am% — Am 2 N (s + s') — uu + tt' + s'u + su + 

(49) 

3ss') + 2 (g z R N ) 2 (m 2 N - s)(m 2 N - s') + 2 (g z - N ) 2 (m 2 N - u)(m 2 N - u') + 



91 (Ag R m\ t + f 2 (0)(Am% - Am 2 N (u + u) + 3uu' + tt' + s'u + su - ss')) 

First part Re(M' arf M' a Z ) l,r of the interference term Eq. (|23|) has the following structure for 
the left-handed incident electrons 
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Re(M'M') L 



8« 3 tt 3 (1 - 2s 2 w ) 



c w s w 



g z L - N (2(l + F 2 (0))m%+ 



2(t - (1 + F 2 (0))(s + s'))m 2 N + s(-s' + 2s' + -u') + -F 2 (0)(tt' + s'u - uu')) + 



g z R N (2(l + F 2 (0))m^ + 2(t - (1 + F 2 (0))(u + u'))m^ + 2™' + -F 2 (0)(us' - ss' + 
tt' + su + 3uu')) + f 2 (0)(Am% - 2(s + s' + u + u')m^ + + ( s + + u ')) + 



(50) 



/ 2 (0)F 2 (0) 



Am 



{16m% - 8(s + s' + u + u') m JV + 4 ( s + + u')m 2 N + i'(tt' - (s - u)(s' - u'))) 



N 



and for right-handed electrons we have 



Re{M an M atZ ) R = 

cw 



(2(1 + F 2 (0))m%+ 



2(t - (1 + F 2 (0))(s + s'))rr& + s(-s' + 2s' + -u') + -F 2 (0)(tt' + s'u - uu')) + 



gt N (2(l + F 2 (0))m% + 2(t - (1 + F 2 (0))(u + u'))m 2 N + 2uu' + -F 2 (0)(us' - ss' + 
tt' + su' + 3uu')) + / 2 (0)(4tt4 - 2{s + s' + u + u')m^ + tt' + (s + u)(s' + u')) + 



(51) 



M^JM (16< 8(s + s' + u + u')< + 4( S + + u')< + - (. - u ){ s' - u'))) 



For simplicity, we have introduced a set of coupling constants g z R N defined as 
9 Z r7l N = /i(0)±5i(0), 



(52) 



1 Z-N 



ie 



g z R N ^w + + gf-fw. + t^^ p (h ~ k 2 ) p / 2 (0) 



2mjv 



where t*7± = are the chirality projector operators. 
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